Polarizibilities as a test of localized 
approximations to the self-interaction 

correction 



J. Messud°'^ * Z. Wang"'^'^ P. M. Dinh"'^ 
P.-G. Reinhard^ and E. Suraud"'^ 

" Universite de Toulouse; UPS; 
Laboratoire de Physique Theorique (IRSAMC); F-31062 Toulouse, France 

^ CNRS; LPT (IRSAMC); F-31062 Toulouse, France 

Institute of Low energy Nuclear Physics, Beijing Normal University, Beijing 

100875, China 

'^Institut fiir Theoretische Physik, Universitdt Erlangen, D-91058 Erlangen, 

Cermany 



Abstract 

We present applications of the recently introduced "Generalized SIC-Slater" scheme 
which provides a simple Self-Interaction Correction approximation in the frame- 
work of the Optimized Effective Potential. We focus on the computation of static 
polarizabilities which are known to constitute stringent tests for Density Functional 
Theory. We apply the new method to model H chains, but also to more realistic 
systems such as C4 (organic) chains, and less symmetrical systems such as a Nas 
(metallic) cluster. Comparison is made with other SIC schemes, especially with the 
standard SIC-Slater one. 
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1 Introduction 



Density-functional theory (DFT) has become over the years one of the most 
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powerful theories for the description of complex electronic systems ranging 
from atoms and molecules, to bulk solids. It allows realistic calculations of 
an ever increasing number of systems in physics and chemistry [T][2][5] . As the 
exact functional is not known, most applications employ the Local Density 
Approximation (LDA), see e.g. |5, or its extension to the Generalized Gradi- 
ent Approximation (GGA) [5]. In spite of their successes, these approaches still 
have deficiencies. In particular, the self-interaction error spoils single-particle 
properties as, e.g., the ionization potential Pf7] . Another critical detail where 
LDA and GGA usually fail is the polarizability in chain molecules [SE]- An 
intuitive and efficient solution is to augment LDA by a Self-Interaction Cor- 
rection (SIC) [TOfTT] . i.e. to introduce an explicit orbital dependence of the 
functional by subtracting by hand the spurious self-interaction. The drawback 
is that it produces a state-dependent mean-field Hamiltonian which requires 
extra efforts to enforce orthogonality of the single particle basis [TT]|12II13|[T^ . 
The optimized effective potential (OEP) method [T^16II17|[TB] overcomes that 
complication as it allows to define the best common (state-independent) local 
mean-field potential V"(r). Indeed some crucial features of the underlying SIC 
as, e.g., the localization of states or the derivative discontinuity p^][20j are 
found to be maintained through an OEP procedure ^Ij, for a comprehensive 
overview see [12] . But the exhaustive SIC-OEP equations are difficult to handle 
and are thus often simplified. A most popular approximation is the so-called 
Krieger-Li-Iafrate (SIC-KLI) approach p!6p!7] and, in a further step of simpli- 
fication, the SIC-Slater approximation [22]. However, SIC-KLI and SIC-Slater 
approximations can easily miss crucial features of SIC as, e.g., the performance 
with respect to polarizability, even if accurate exchange-correlation potentials 
are used [23] (note that the same conclusion holds at the exact-exchange level 

A key question is the localization of orbitals, as had been already observed 
rather early [2l], which then minimizes the electronic interaction energy. In- 
deed states which optimize SIC tend to be localized whereas states emerging 
from a common local mean-field tend to more delocalization. These contradict- 
ing demands can be bridged by full OEP but spoil approximate treatments 
as SIC-KLI or SIC-Slater approximation. A fairly convenient way out is to 
use two sets of orbitals. That was already proposed in [2S] and has been used 
to improve on the SIC-KLI approximation [25]|27f28] . although not yet in a 
fully consistent and variational form. Taking up the double-set idea, we have 
recently proposed a SIC-OEP scheme relying on two sets of complementing 
orbitals [29]. At the purely stationary full SIC level, it is mostly a matter of 
convenience, while a double-set technique becomes crucial in practical appli- 
cations of the time-dependent SIC [T4] . In this paper, we will also demonstrate 
that this method can be powerfully exploited in stationary calculations in the 
frame of SIC-OEP. Indeed the two sets are connected by a unitary transfor- 
mation, thus building the same total density. One of the sets remains spatially 
localized while the other set is free to accommodate a common local potential 
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and/or minimal energy variance. The spatially localized set is determined by 
variation of the SIC energy with respect to unitary transformation coefficients 
which allows that these states fulfill what is called the "symmetry condition" 
P^25f29ll3U] . The localized set shows a much better performance with respect 
to standard SIC-KLI or SIC-Slater approximation. The resulting formalism is 
called "Generalized SIC-OEP" [29]. We further simplified the resulting equa- 
tions, following the track of the SIC-Slater approximation and developed in a 
strictly variational manner a double-set treatment of the SIC-Slater approxi- 
mation to OEP, namely the "Generalized SIC-Slater" approximation. It is to 
be noted that a very similar development is found in [21]. As one of the sets 
remains spatially localized, it validates the Generalized SIC-Slater approxi- 
mation to Generalized SIC-OEP built from this set, while maintaining key 
features of the full SIC scheme: it is energetically advantageous for the SIC 
energy, permits to re-establish potential energy surfaces (PES), and performs 
fairly well in the calculations of polarizabilities. We pointed out that the Gen- 
eralized SIC-Slater emerges naturally because of the localization of one set 
of orbitals [29], whereas Korzdorfer et al. [31] compared Generalized SIC-KLI 
and full Generalized SIC-OEP with standard SIC-KLI and OEP and also with 
the approximate Generalized SIC-KLI of [28j . It is to be noted that the prac- 
tical introduction of localized orbitals within the SIC-KLI approximation has 
also been proposed in [26ll27f28j without explicit "symmetry condition" . The 
form of the SIC-KLI approximation used in those papers is not exactly what 
comes out from the more fundamental "Generalized SIC-OEP" [31]. After a 
brief presentation of the formalism, we apply it to model hydrogen chains, and 
to various other systems such as C (organic) chains and Na (metallic) clusters. 



2 The Generalized SIC-Slater formalism 

2. 1 Summary of SIC equations 

We briefiy summarize the formalism using for simplicity a notation without 
explicit spin densities. The generalization to these is obvious. The calculations 
later on use, of course, the full spin density functional. 

The starting point for the formulation of SIC is the SIC energy functional for 
electrons 

N N 

EsiC = ^kin + ^ion + ^LDA[p]-Xl ^LDa[P/3] , P = ^ P(3 , Pf} = 1^/3^ (1) 

P=l (3=1 

where E'ldaIp] is a standard LDA energy- density functional (which contains 
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both the Hartree and the exchange-correlation energies in our notations) com- 
plementing the kinetic energy i?kin and the interaction energy with the ionic 
background i^ion- The last term is the SIC correction. The densities and p 
are defined from the set of occupied single-particle states {ipp, [3 = l...A^}. The 
SIC equations are obtained by standard variational techniques within impos- 
ing explicit orthonormalization of the orbitals by a set of Lagrange multipliers 
Aa/3. We now introduce a second set of orbitals {(pi} related to the previous 
one by a unitary transformation within the set of occupied states (i.e. leading 
to the same total density p : ipi = J^aKai^a) which diagonalizes the Xa/s- We 
can then recast the resulting equations in eigenvalue equations 



hsiclfi) = (2) 
= iijf3\Up-UM (3) 



with the SIC Hamiltonian reading 



hsic = /iLDA - 51 I ^« ) I (4a) 

a 

/iLDA = |^ + t/LDA[p](r), With f/LDA [p] (r) = (4b) 

^ 5^^^^^^ r (4c) 

where /ilda is the standard LDA Hamiltonian. The coefficients Uia of the 
unitary transformation for given diagonal orbitals (pi are determined such 
that the ipa satisfy the symmetry condition ([3]). The ipa are called localized 
orbitals because they are spatially much more localized [25f30] . 



2.2 SIC and OEP 



The eigenvalue equation ([2]) employs a non-local Hamiltonian /isic, see Eq. ( I4all . 
which complicates the numerical handling. In [22], we proposed to apply the 
OEP formalism to this two-sets SIC formulation, to find the best local ap- 
proximation to its Hamiltonian. A very similar development is found in |31j . 
We start from a set ipi which satisfies the eigenvalue equations (this set is 
not exactly the same as that of the exact SIC equation because additional 
restriction of the Hilbert space is imposed here - this point being clarified, we 
will employ the same symbol to simplify the notations) : 



/iLDA(r) - Vq{y) v9i(r) = eiipi{Y) , (5) 
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where Vq is a local and state-independent potential which needs to be opti- 
mized to minimize the SIC energy ([1]). It is important to note that this energy 
is still expressed in terms of the ipai linked by a unitary transformation to the 
(pi and which satisfy the symmetry condition ([3]) in our case. The optimized 
effective potential Vo(r) is found by variation 5 Eq,ic / 5Vq{y) = 0. We obtain 
Vo = Vs + Vk + Vc, with [29II3T] 



Vs = Y.^-^U^BAm'] , (6a) 



VK = -T.(Y.\^^\'^>^p)(^^^\yO-Ul^DA[m']m , (6b) 



p 

P.(r) = -;^E^- /dr'(K.(r')-f/LDA[|^a|1(rO)c(r')G.(r,r'),(6d) 
(r) a V / 



2.3 Generalized SIC-Slater 



The involved OEP equations can be simplified by exploiting the property that 
the ipa are spatially localized |25], which yields Vb|'?/'a) ~ t^LDA[|'^aP]|V'a)- This 
allows to employ the SIC-Slater approximation to OEP, yielding [29] : 



^o(r)-E^^^LDA[|^.r](r) , (7) 

Note that this equation has the form of a SIC-Slater approximation [2T][32] 
but is constructed from the locahzed orbitals ifja and is applied to the diag- 
onal orbitals ipi. We called this new scheme "Generalized SIC-Slater" (GSlat) 
approximation, which differs from the standard SIC-Slater scheme because of 
the two basis sets involved here and which, therefore, has more flexibility. The 
practical scheme for GSlat can be summarized as follows : i) Eq. ([5]) gener- 
ates the "diagonal" set ipi of occupied states; ii) the unitary transformation 
serves to accommodate the symmetry condition ([3]) which, in turn, defines the 
"locahzed" set ipa] Hi) the latter set enters the OEP Vq as given in Eq. ([71). 



2.4 Generalized SIG-KLI and approximations thereof 
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Even if the GSlat approximation might be satisfying in most cases, it is worth 
discussing in more detail the SIC-KLI correction (16bp . The use of the local- 
ization of the ipa allows to keep only the diagonal terms a = f3, that is, 

^ - E (E K^^I')(^"I^0 - U^DAlMl^Pa). (8) 
P a ^ i ' 

In contrast to the GSlat approximation, fluctuations of VolV'a) around t/LDA[|V^a 
are not neglected, even if one expects them to remain small. 

It is now instructive to use the following identity : 

I E^raV'iP = E I^LV'iP + EKa<^i)Ka<^i)% 

% % i,j 

to multiply it by (V'alVo — f/LDA[|V'aP]|V'a) and to sum over a. We thus obtain : 
where 

Aij = ipiip*'^V*^Vjaii^a\Vo - t/LDA[|V'«ni^«)- (10) 

a 

First, we have the expression Vk ~ ^ Z^i where we used the approximate 
form ([8]) of the SIC-KLI correction. If one further assumes that for i j, 
\Aij\ ^ l^iil, one obtains the following approximation of the SIC-KLI correc- 
tion 

VK--T.\^o.\\^a\Vo-Ui^BA[m'Ua) , (H) 
P a 

which is the form proposed in [26|27|28] without reference to the symmetry 
condition. This expression can thus be derived from an approximation of the 
so-called "Generalized SIC-KLI" potential. We will call it "Locahzed SIC-KLI" 
(Loc. SIC-KLI) thereafter. 

The justification of the approximation \Aij\ <^ l^jjl remains to be clarified. 
A way to circumvent the formal difficulty is to consider practical applications 
to see how the approximation performs. Still, the assumption probably holds 
for spatially symmetric systems, as will be discussed lengthly when tested 
on hydrogen chains (see Sec. 13.11) . The case of asymmetrical systems however 
remains to be explored in more detail, especially in comparison to the formally 
better founded forms (l6b1) or ([H]). 
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2. 5 Summary 



All SIC mean-field Hamiltonians presented above enter a Schrodinger-like 
equation of the form h\ipi) = ei\ipi). We have thus summarized them in table 
[H so that one can easily track the various contributions from one Hamiltonian 
to another. Note that the symmetry condition ([3]) should be added for the last 



Tab: 



Expression of h in h\ipi) = ei\(pi) 


Method 


/ilda[p] 


LDA 


hhDAip] - f^LDA ^ 




Average Density SIC 


j ^ 


Standard SIC-Slater 


/iLDA[p] ^ ' ?7LDA['0a'] 
a ^ 


Generalized SIC-Slater 


hLDA[p] '"^"'VldA [IV'aP] 
a ^ 

--y2\^a\\^a\Vo-ULDA[\i^a\''Ma) 
Pa 


Localized SIC-KLI 


hhDAip] - I]q f^LDA 




Exact SIC (benchmark) 


le 1 



The hierarchy of mean-field Hamiltonians, from simple-most LDA (top line) to full 
SIC (bottom line), called "exact" SIC. 



three schemes, to define the localized states i/ja required in the corresponding 
Hamiltonians. 

2.6 Computational details 

The test cases presented in this paper have been obtained using a full 3D DFT 
code originally developed for large scale calculations of the dynamics of metal 
clusters [33] and later on smah hydrogen clusters [3l], and now extended to 
treat any organic system [35]. The electronic wave functions are represented 
on an equidistant 3D grid with fixed mesh size (between 0.4 and 0.8 Cq, de- 
pending on the studied system). The ionic background is treated by means of 
pseudopotentials, using either local ones (Na, H) [36] or non-local ones to treat 
organic systems [27] • Even in the case of hydrogen, we use a pseudopotential 
in order to regularize the Coulomb singularity at origin (the Giannozzi pseu- 
dopotential as in [23]). The pseudopotential parameters, especially the core 
size, fixes the optimal grid representation. For the molecular chains calcula- 
tions, we use boxes of 40 — 52 qq in the longitudinal direction (according to the 
case) and 20 Qq in the transverse directions. For the Na5 cluster, we use a box 
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of 38 Oo in each direction. Those sizes are sufficiently large to provide good 
convergence properties of calculations. For completeness, computations have 
been checked using larger boxes with no noticeable differences in the results. 
Electronic ground states are obtained by damped gradient iterations. 



3 Static polarizability results 

We had shown in [29] that GSlat solves the problem with potential-energy 
surfaces encountered in the standard SIC-Slater scheme and produces good 
results for the polarizability of the C atom. Here we will show through full 3D 
calculations that it also yields favorable results for more complex structures: 
model H chains, C4 chains and a Nas cluster. We compare the GSlat results 
to LDA, ADSIC (Average Density SIC) [38], standard SIC-Slater and exact 
SIC results, the latter being the benchmark. For the comparison, we use the 
static dipole polarizability as a most sensitive test for DFT approximations 
[39] . Considering a system put inside an electrical field E, the polarizability is 
defined as Ui = dfii/dEi, where /ij is the dipole moment along the i direction 
and Ei the electric field along i. 



3.1 Hydrogen chains 

Linear chains of H atoms constitute (highly) simplified model systems for 
various important chain or chain-like molecules such as in particular poly- 
acetylene with its remarkable properties [ID] . These model systems are of great 
interest to investigate DFT schemes [2S|^ as they are particularly difficult 
to be correctly described within LDA [27]. They thus provide a critical test. 

Our calculations on polarizabilities in hydrogen chains are presented in Fig. [1] 
and are compared with previous results [28ll23f42j . For sake of a fair compari- 
son, we have taken care of using the same pseudopotential as used in former 
calculations because we found that using different pseudopotentials leads to 
different absolute values of the highly sensitive polarizabilities. This point 
may look surprising especially in the case of hydrogen for which inserting a 
pseudopotential is a matter of practical convenience for regularizing grid rep- 
resentations. However, as the grid size is optimized to the regularizing core, 
different core widths lead to slightly different finite representations of the wave 
functions (having of course the same energy) . This may deliver slighly different 
values of polarizabilities. This has to be kept in mind when comparing with 
MP4 results which are computed in a different fashion. But the comparison 
between the various DFT approaches stays on safe grounds. 
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4 6 8 10 12 

Chain length n 



Fig. 1. Top : Longitudinal polarizabilities of H„ chains (per hydrogen atom, in 
ao^) as a function of length n, for an alternation of 2 and 3 ao bond length in 
various calculations; we present our own LDA, SIC, standard SIC-Slater and GSlat 
results (open symbols) in comparison with calculations from other groups: standard 
SIC-KLI [23], standard SIC-OEP [23], Loc. SIC-KLI [28] and quantum chemistry 
MP4 [42]. Bottom: Transverse polarizabilities of H chains (our own results only). 



We first mention that our LDA and SIC calculations perfectly match pre- 
viously published results [13] (not shown in the figure). Furthermore the 
trends are rather systematic: LDA strongly overestimates polarizabilities, as 
expected, and SIC comes much closer to MP4 results. Still the most relevant 
comparison, in our opinion, is that between SIC and approximations thereof, 
because of the pseudopotential effects mentioned above. 

There are several interesting points showing up from this comparison. First, 
while standard SIC-Slater (open circles) and standard SIC-KLI (full circles) 
calculations lead to results of varying quality, we see that GSlat (open squares), 
even if not perfectly matching exact SIC, leads to overall acceptable results 
showing more regular and realistic trends. However the agreement in absolute 
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values tends to degrade with increasing chain length, in relative values (of 
GSlat compared to SIC) : from +8% to +25% from H4 to H12. Hence one 
can wonder whether less dramatic approximations to the Generalized SIC- 
OEP formalism would improve the results. An obvious next step is to use 
the Generalized SIC-KLI correction instead of GSlat. We show in Fig. [T]the 
Loc. SIC-KLI results obtained by [28j with the addition of the approximate 
Generalized SIC-KLI term as given in Eq. (11 II) . Mind that these results used 
another localization criterion than the symmetry condition. As seen in the 
figure, perfect agreement with SIC is achieved (compare large open triangles 
with black diamonds). This thus validates a posteriori the assumption done 
to obtain (fTTl) . We however recall that Eq. (fTTj) has no robust fundation. The 
Loc. SIC-KLI results nevertheless demonstrate the promising possibility to use 
a numerically less costly localization criterion than the symmetry condition. 
Indeed, the equivalent GSlat in ^28j perfectly agrees with ours. 

The mechanism invoked to explain the better performance of SIC-KLI (stan- 
dard and generalized) on polarizabilities, over mere LDA or standard SIC- 
Slater, is that the SIC-KLI correction, or the so-called response part of the 
exchange-correlation potential defined in [S] , produces a counter-field effect [8] : 
with an external field, the exchange-correlation potential behaves globally 
against it. Electronic motion is hindered and polarizabilities are thus reduced 
compared with a LDA treatment. The counter- field effect is all the more effi- 
cient when one goes from SIC-KLI to OEP [^TpS] or in a Generalized SIC-KLI 
treatment [28j. In GSlat, no response potential is present, so no counter-field 
effect is expected here. Indeed, the average GSlat exchange-correlation poten- 
tial in the presence of an external field is not opposite to the electric potential, 
as the SIC-KLI or the OEP ones are (see e.g. bottom panels of Fig. 2 in [28] 
which correspond to our GSlat but obtained with another localization criterion 
than the symmetry condition). Actually, the performance of our GSlat approx- 
imation stands in the property that much higher barriers (than in standard 
SIC-Slater or SIC-KLI) between H2 units appear in the exchange-correlation 
potential, which hinder more the electronic motion. This has to be linked to 
the localized character of the orbitals that enter into its calculation |28j. 

It would be extremely interesting to also test the Generalized SIC-KLI and 
its various approximate forms in less symmetrical systems. This calls for a 
systematic study which will be reported in a forthcoming paper. Keeping 
this in mind, we will nevertheless present in the following a few examples of 
applications to less symmetrical systems in the case of GSlat in comparison 
to exact SIC. 

Before doing so, let us remark that the case of full standard OEP results 
of [23] (close squares) deserve a special comment. Indeed, up to fiuctuations, 
the obtained results perfectly match the MP4 ones (stars), and thus somewhat 
differ from SIC ones. There is here a matter of interpretation in the sense that. 
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if the aim is to match MP4 results, the results are perfect. But if it is to find 
a good approximation to SIC then the agreement is not dramatically better 
than GSlat ones (relative values of OEP compared to SIC : —6% to —13% 
from H4 to H12) and is worse than the Loc. SIC-KLI ones. Moreover GSlat 
and Loc. SIC-KLI calculations seen as approximations to full Generalized SIC- 
OEP are cheaper schemes, easily applicable to much more complex systems 
than hydrogen chains. 

Note finally that Fig. [1] also shows the transverse polarizability for those hy- 
drogen chains. In that case, the GSlat approximation reproduces perfectly the 
exact SIC results, as expected. 

Now that the capability of GSlat on hydrogen chains has been checked, we 
deeper analyze the properties of the system in a small chain, namely II4. We 
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Distance between the center of mass of each Hg (ag) 



Fig. 2. Longitudinal (top) and tranverse (bottom) polarizabilities of H4 chains, 
according to the H2-H2 center of mass distance, for various SIC schemes as indicated. 



present in Fig. [2] the values of the (longitudinal and transverse) polarizabili- 
ties of H4 chains, according to various H2-H2 center of mass distances. Here 
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we used the experimental value of the H2 bond length, that is 1.46 oq. The 
data labeled "SIC" constitutes our benchmark. LDA (stars) overestimates po- 
larizabilities which was expected on the ground that LDA has a tendency to 
overmetallize bonding. The simplified ADSIC (open squares) scheme [3H] gives 
in general rather poor results. Mind that ADSIC nevertheless allows a fair re- 
production of bonding properties in poly-acetylene [H]. It obviously fails in 
the case of the more sensitive polarizability. GSlat (crosses) in turn repro- 
duces very well the exact SIC tendencies, while the standard SIC-Slater (full 
squares) is completely wrong for intermediate intermolecular distances. This 
mismatch is correlated to a similar failure of standard SIC-Slater in the po- 
tential energy surface at intermediate distances [29] . And both failures can be 
tracked back to delocalization effects of the orbitals at critical configurations. 
GSlat allows to keep the wave functions entering the Hamiltonian localized and 
thus performs much better. We checked that standard SIC-KLI (not shown 
in the figure) does not cure this mismatch. For large intermolecular distances, 
standard SIC-Slater and GSlat results come close to each other because there 
remain two separated H2 molecules, which have each only one electron in each 
spin subspace. 



3.2 The C2 dinner and the C4 chain 

Another interesting case is provided by small carbon chains whose electrical 
excitation properties are well studied [1^461147] . We consider here two ex- 
amples, namely the C2 dimer and the C4 chain. Since, to the best of our 
knowledge, there exist no experimental results for the polarizabilities of those 
systems, our aim is to compare various theoretical approaches using exact SIC 
as a benchmark. We recall that we already demonstrated the quality of GS- 
lat in the case of a single carbon atom in [22] for exchange only calculations. 
In exchange correlations calculations, the C atom polarizability is again well 
reproduced by GSlat. As the electronic cloud in the C atom is slightly axially 
deformed because of the single occupation of and 2py orbitals (while the 
2pz orbital is unoccupied), one measures different polarizabilities along the 
X, y axes and along the z axis. To put the subsequent results on C molecules 
into perspective, we quote here briefly the polarizations for the C atom : along 
z axis, az = 10.40 Gq^ for both SIC and GSlat, along x,y axes, a^^y = 11.52 
ao^ for GSlat and 11.76 for SIC. 

In Fig. [3l we present the longitudinal and transverse polarizabilities for C2 
and C4 calculated in various approximations. The calculations of [ISMH] yield 
comparable values. In [15], the longitudinal polarizability for C2 is ay = 25 Oq^ 
for the ab initio methods and 34 a^^ for LDA/GGA, while the transverse one is 
a± =25 Qq^ or 100 respectively, the latter value being a strange exception. 
The results for C4 are a\\ =92 or 94 oq^ and a± = 30 or 32 ao^. Our results are 
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Fig. 3. Transverse and longitudinal polarizabilities of the C2 molecule (left) and the 
C4 chain (right), calculated in various SIC schemes. Horizontal lines emphasize the 
SIC benchmark values and ease the comparison with the other results. 

generally lower for aj^. It is worth noting that our calculations differ in the 
employed functionals and pseudopotentials which both can have a sensitive 
influence on the results. Thus the comparison as a whole looks satisfying. To 
stay on the safe side, we concentrate on the comparison of approaches within 
the same setup. 

Fig. [3] shows again that GSlat provides a very good approximation to exact 
SIC in C2 and much better than any other approximation. The situation is 
more mixed in the case of C4. The values are larger and the relative effects are 
smaller than for the C2 dimer. GSlat comes still closest to SIC for the longitu- 
dinal mode, but standard SIC-Slater is competitive for the transverse mode. 
Still, when considering all cases together (C2 and C4, transverse and longitudi- 
nal polarizabilities), it is clear that GSlat provides a very good approximation 
(generally the best one) to the exact SIC. 

3.3 Metal clusters 

As a final test case, we consider a small sodium cluster representative of simple 
metallic systems. We have chosen the Nas cluster because it has a very soft 
electron cloud and is thus a most critical test case amongst metallic particles 
[l9] . The cluster is planar (see structure in the insert of Fig. H]) which corre- 
sponds to a triaxial shape, and has accordingly different polarizabilities along 
the three major axes of the system. Fig. H] shows the polarizabilities of the 
Nas. We obtain much larger absolute values of polarizabilities than in the case 
of organic systems due to the metallic nature of bonding (delocalization and 
lower binding). Not surprisingly, LDA performs rather well, for sure better 
than in organic systems, as is to be expected for a simple metallic system. 
Even if LDA works very well on those kind of systems, we still see a non neg- 
ligible difference with the benchmark SIC for the X direction, about 7% for 
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Fig. 4. Polarizabilities of the Nas metal cluster (displayed in the right panel), for 
various SIC schemes as indicated. Horizontal lines emphasize the SIC benchmark 
values and ease the comparison with the other results. 



LDA. Again, GSlat stays closest to the benchmark in all cases. 



4 Conclusion 



We have tested a newly developed DFT-SIC scheme, called Generalized SIC- 
Slater (GSlat), with respect to polarizability in chain molecules and a soft 
metal cluster. GSlat starts from the Optimized Effective Potential (OEP) ap- 
proach to SIC and handles that in terms of two different sets of single- 
particle wave functions. One set is taken for the solution of the OEP equa- 
tions, thus diagonal in energy and most likely delocalized. The other set is 
used in setting up the SIC energy which becomes lowest for localized wave 
functions. Both sets are connected by a unitary transformation which leaves 
key features as, e.g, the total density invariant. Using that double set allows 
to accommodate two conflicting demands, energy diagonality versus locality. 
The unitary transformation is determined by minimization of the SIC energy 
which leads to what is called the symmetry condition, a key building block of 
the SIC equations. The localized character of the SIC optimizing set is well 
suited to justify the steps from OEP to SIC-KLI and further to the SIC-Slater 
approximation. Thus SIC-OEP with double-set representation and subsequent 
SIC-Slater approximation leads to the GSlat scheme. By virtue of the double- 
set technique, it has more flexibility than standard SIC-KLI or SIC-Slater 
approximation. 

As it is known that the polarizability in chain molecules is a sensitive observ- 
able for DFT approaches, we have investigated the performance of the new 
scheme with respect to polarizability in a variety of critical test cases : H 
chains which mimic the electronic properties of polymers, the C2 dimer, a C4 
chain, and Na5 as a soft small metallic particle. The results demonstrate that 
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the GSlat approximation comes generally close to the values from exact SIC 
which we use here as a benchmark. It solves the pathologies of the standard 
SlC-Slater approximation which occur in critical (transitional) molecular con- 
figurations such that its performances depend much less on the kind of studied 
system and configuration (which is not the case in standard SIC-Slater or SIC- 
KLl). For long H chains, some deviation is observed, which remains reasonable 
compared to the other schemes and to the less important numerical cost of 
the GSlat scheme. However addition of the SlC-KLl correction (coming from 
Generalized SIC-OEP) allows to improve substancially the results. Neverthe- 
less, no strong argument justifies a priori the approximate form in Eq. flTT]) . 
even for spatially symmetrical systems. And the case has yet to be tested for 
asymmetrical systems. We here checked that GSlat performs fairly good as 
well in H chains than in less symmetrical systems as C chains or Na clusters. 
We finally checked that the improvement of the SlC-KLl correction over GSlat 
is negligible when looking for other observables as energies. 

This work was supported, by Agence Nationale de la Recherche (ANR-06- 
BLAN-0319-02), the Deutsche Forschungsgemeinschaft (RE 322/10-1), and 
the Humboldt foundation. 
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